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The cosmic acceleration is one of the most significant cosmological discoveries over the last century. 
The two categories of explanation are exotic component (dark energy) and modified gravity. We 
i— ( ' constrain the two types of model by a joint analysis with perturbation growth and direct H(z) data. 

Though the minimal \ 2 of the ACDM is almost the same as that of DGP, in the sense of consistency 
we find that the dark energy (ACDM) model is more favored through a detailed comparison with 
' the corresponding parameters fitted by expansion data. 
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The acceleration of the universe is one of the most significant cosmological discoveries over the last century 
Various explanations of this acceleration have been proposed, see Q for recent reviews with fairly complete lists 
l— ~~ '■ of references of different models. However, although fundamental for our understanding of the universe, its nature 
| remains as a completely open question nowadays. 

There are two main categories of proposals. One is that the acceleration is driven by some exotic matter with 
' negative pressure, called dark energy. The other suggests that general relativity fails in the present Hubble scale. 

The extra geometric effect is responsible for the acceleration. Surely, there are some proposals which mix the two 
1 categories. Mathematically, in the dark energy model we present corrections to the right hand side of Einstein equation 
(matter part), while the correction terms appear in the left hand side of Einstein equation (geometric part). 
\Q • ACDM model is the most popular and far simple dark energy model, in which vacuum energy with the equation 
of state (EOS) w = — 1 accelerates the universe. From theoretical considerations and by observational implications, 
people put forward several other candidates for dark energy, such as quintessence (— 1 < w < —1/3), phantom 
(w < —1), etc. Also there are many possible corrections to the geometric part of the theory. One of the leading 
modified gravity model is Dvali-Gabadadze-Porrati (DGP) model Q, for a review, see In the DGP model, the 
bulk is a flat Minkowski spacetime, but an induced gravity term appears on a tensionless brane. In this model, gravity 
appears 4-dimensional at short distances. But, at a distance larger compared to some freely adjustable crossover scale 
ro it is altered through the slow evaporation of the graviton off our 4-dimensional brane world universe into an unseen, 
yet large, fifth dimension. 

We should find the correct, at least exclude the incorrect models in the model sea. The first step is to discriminate 
between dark energy and modified gravity, whose nature are completely different. To construct a model simulating the 
accelerated expansion is not very difficult. That is the reason why we have so many different models. Recently, some 
suggestions are presented that growth function S(z) = Sp m / p m of the linear matter density contrast as a function of 
redshift z can be a probe to discriminate between dark energy and modified gravity models. The growth function 
can break the degenerations between dark energy and modified gravity models which share the same expansion history. 

There is an approximate relation between the growth function and the partition of dust matter in standard general 
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1.46 ±0.29 
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TABLE I: Observed perturbation growth as a function of redshift z, see also fl5| 

relativity 0, 

/-^ = n^. CD 

a ma 

where f2 m is the density partition of dust matter, a denotes the scale factor, and 7 is the growth index. This relation 
is a perfect approximation at high redshift region. Also, it can be used in low redshift region, see for example [8J. It 
is shown that the relation ([T]) is also valid in the case of modified gravity theory [j| . The theoretical value of 7 for 
ACDM model is 6/11 @, ||, while, for spatially flat DGP model is 11/16 0. The observation data of perturbation 
growth are listed in Table 1. 

It is shown that ACDM model is consistent with the current growth data [l(| • The data seem to weaken a spatially 
flat DGP model, whose 7 = 11/16. However, it is found that the growth index 7 is 4/7 in a non spatially flat DGP 
model [ll|, which is very closed to the index of ACDM. Thus, the DGP model may be still consistent with current 
growth data. 

In this article we take a different strategy. The previous works were concentrated on the limit of the growth index 
and made some approximations on it (often the high z limit was assumed and an approximation was made at linear 
order), in which only approximate asymptotic value of 7 can be obtained. In fact, the perturbation growth / is a 
variable with respect to z, as displayed in Table 1. By using these growth data we constrain the parameters in ACDM 
model and DGP model, respectively. 

The other one which is very useful but not widely used in model constraint data is the set of direct H(z). H(z) 
is derived by a newly developed scheme to obtain the Hubble parameter directly at different redshift [17| . which 
is based on a method to estimate the differential ages of the oldest galaxies [l8j]. By using the previously released 
data [l9|, Simon et al. obtained a sample of direct H{z) data in the interval z € (0,1.8) [20], just as the same 
interval of the data of luminosity distances from supernovae. For the present sample of growth data derived with the 
assumption of the expansion behaviors of the universe, we will present joint fittings to obtain the constraints on the 
ACDM and DGP, respectively. Then, through comparing with allowed regions by expansion constraint using different 
observations, including supernovae (SN), cosmic microwave background (CMB), baryon acoustic oscillations (BAO) 
etc., we examine which model is more self-consistent. 

This article is organized as follow: In the next section we construct the evolution equation for / in a very general 
frame. In section III, by using the growth data and H{z) data we present the parameter constraints of ACDM and 
DGP, respectively. Our conclusion and some discussions appear in the last section. 

II. THE EVOLUTION EQUATION FOR THE GROWTH FUNCTION / 

We consider a mixed model in which dark energy drives the universe to accelerate in frame of modified gravity. For 
FRW universe in modified gravity, the Friedmann equation can be written as, 

7-7-2 ^ / ■ ~s 87rG 

H + — + h(a,a,a) = ——{p m + p e ), (2) 
a z 6 

where H denotes the Hubble parameter, h comes from the corrections to general relativity. p m and p e represent the 
density of dust matter and the exotic matter, respectively. A dot implies the derivative with respect to cosmic time 
t. Comparing with the corresponding Friedmann equation in standard general relativity, we obtain the density of 
effective dark energy, 
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Here we call h geometric sector of dark energy. The behavior of the effective dark energy has been separately discussed 
in some previous works. For example, it is investigated in detail in a modified gravity model where a four dimensional 
curvature scalar on the brane and a five dimensional Gauss-Bonnet term in the bulk are present [2l[ . 

For any modified gravity theory, Bianchi identity is a fundamental requirement. Using the continuity equation of 
the dust matter and the Bianchi identity, we derive, 

p de + 3Hp de (l+w de ) =0, (4) 

which yields, 

ldlnpde feL s 

w de = -1 - ^-Tj . (5) 

3 a ma 

where Wde is the EOS (equation of state) of the effective dark energy. 

After the matter decoupling from radiation, for a region well inside a Hubble radius, the perturbation growth 
satisfies the following equation in standard general relativity [22| . 

8 + 2HS - iTvGp m 8 = 0. (6) 

It is found that the perturbation equation is still valid in a modified gravity theory if we replace the Newton constant 
G with an effective gravitational parameter G e f t , which is defined by Cavendish- type experiment |23| (This point 
may need more studies.). G e // may be time-dependent, for example in Bran-Dicke theory and in generalized DGP 
theory [Hj]. 



With the partition functions, 
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the perturbation equation ^ becomes, 



(In 6)" + (\n6Y 2 +^2 + ^pj (In 6)' = ~aQ m , (10) 

where a prime denotes the derivative with respect to In a, a is the strength of the gravitational field scaled by that of 
standard general relativity, 



rim and Qk redshift as 



n = ^°( 1 + ^) 3 (12) 

m n m0 (i + ^)3 + n fe0 (i + 2 )2 + 8^' 

" fi m0 (l + 2 )3 + ^ fc0 (l + z )2 + 8^' 



where denotes the present value of a quantity. The growth function defined in (TTJ is just (In 5)' . Thus (fTQ|) generates 

/' + f + 



1 3 

-(l + n k ) + -w de (n m + n k - 1] 



f = \atl m , (14) 



where we have used 



— = -Q k - - [fl m + (1 + W d e)ttde] , (15) 
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and 



Q m + fl k + Qde — 1< 

In ACDM model, we have Wd e = — 1 and a = 1. For the self-accelerating branch of DGP model [Hj ]. 

-1 + fi* 



i + fi m — n k 



and 
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r c is another important parameter in DGP model , which is defined by the relative strength of five dimensional gravity 
to four dimensional gravity r c = G5/G. Here G5 is the five dimensional gravity constant. We define the partition of 

7> as 



One can derive the following relation from Friemann equation of DGP model, 

, n 2 



1 
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(19) 



(20) 



With (Ti"4f and the observed data of / in Table 1, we can fit parameters of the models, either dark energy or modified 
gravity. 



III. JOINT ANALYSIS WITH THE GROWTH DATA AND THE DIRECT H{z) DATA 

In this section we fit f2 m o, f^o in ACDM model and Q m o, ^r c in DGP model with the growth data and direct H(z) 
data by % 2 statistics, respectively. Before fitting with the two sets of data, we present some discussions about them. 

The present growth data in Table 1 are far from being precise. We have a sample consisting of only six points, 
and the error bars of the growth data are at the same order of the growth data themselves. The reason roots in the 
method by which we derive the data set. 

In the present stage we do not find any absolute probes to the perturbation amplitude. People extract the infor- 
mation of perturbation growth from galaxy clustering data through redshift distortion parameter j3 observed in the 
anisotropic pattern of galactic redshifts. We need the galaxy bias factor b to get the perturbation growth / = b(3. 
The current available galaxy bias can be obtained mainly in two ways. The most popular method is to refer to the 
simulation results of galaxy formations 0, US (3 ■ At. the present stage the simulations we obtained only in frame 
of ACDM model. The second method to get the galaxy bias depends on the CMB normalization [ll[. Also, ACDM 
model is involved. Further, to convert from redshift z to comoving distance one should assume a clear relation between 
distance and redshift. For instance Tegmark et al. adopt a flat ACDM model in which f2 m Q = 0.25. They also 
tested that if a different cosmological model is assumed for the conversion from redshift to comoving distance, the 
measured dimensionless power spectrum is varied very slightly (<1%) [25| . 

Hence, we see that people obtain data in Table 1 always by assuming a ACDM model. Its reliability may decrease 
when we use it in the scenarios of other models. Fortunately, it is pointed out that this problem can be evaded at 
least in the DGP model since the expansion history in DGP with proper parameters is very similar to that of ACDM 
[l5l |. Since the the growth data are derived with some assumptions of expansion history, we should fit the model by 
growth data together with expansion data. 

The direct H(z) data are independent of the data of luminosity distances and reveal some fine structures of H(z). 
They have not been widely used in the constraints on dark energy models up to now. Here we present a joint fitting 
of ACDM and DGP with perturbation growth data and direct H(z) data. 

We show the sample of H(z) data in table II. 

Table II displays an unexpected feature of H(z): It decreases with respect to the redshift z at z ~ 0.3 and z ~ 1.5, 
which is difficult to be found in the data of supernovae since the wiggles will be integrated in the data of luminosity 
distances. A study shows that the model whose Hubble parameter is directly endowed with oscillating ansatz by 
parameterizations fits the data much better than those of LCDM, IntLCDM, XCDM, IntXCDM, VecDE, IntVecDE 
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0.09 


0.17 


0.27 0.40 


0.88 


1.30 
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1.53 


1.75 


H(z) (km s^ 1 MpcT 1 ) 


69 


83 


70 87 


117 


168 


177 


140 
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68.3% confidence interval 


±12 


±8.3 


±14 ±17.4 


±23.4 


±13.4 


±14.2 


±14 


±40.4 



TABLE II: The direct observation data of H(z) [H . 



[261 ] . A physical model, in which the phantom field with natural potential ,i.e, the potential of a pseudo Nambu- 
Goldstone Boson (PNGB) plays the role of dark energy, is investigated in [27| • The oscillating behavior of H appears 
naturally in this model. 

For ACDM, in the joint analysis with a marginalization of Hq, \ 2 reads 



'(^mO, Oo) — 



fobs(zi) — fth(Zi] QmO, Oo) 



9 r 



E 



H bs{Zi) — H t h(Zi] fl,nO, Ofeo) 



-i 2 



Hp -72 
0.08 



(21) 



where fobs denotes the observation value of the growth index, and fth represents its theoretical value. We read fobs(zi), 
<Jf ob3 from Table 1 and calculate fth(zi] Q m o, Oo) using (fT4"l) . To get the theoretical value of / using (fT4"]) . we need 
its initial value. Our considerations are as follows. In any dark energy model the universe should behave as the same 
one in some high redshift region such as z — 1000, that is, it behaves as standard cold dark matter (SCDM) model, 
which has been sufficiently tested by observations. In SCDM model we obtain / = 1 by using |T]). So we just take 
/ = 1 as the initial value at high enough redshift region. And the theoretical Hubble parameter reads, 



H th 



Hi [n ro0 (i ± zf + n M (i ± zf + 1 - n m0 - n 



feO 



(22) 



We take the value of present Hubble parameter Hq from the HST key project Hq = 0.72 ± 0.08kms _1 Mpc _1 [3C| . 
The result is shown in fig. 1. 

In DGP model, traditionally, we often use O r< . rather than Oo m fittings. There is no essential difference since 
they are constrained by (|20|) . In the joint analysis with a marginalization of Hq, \ 2 becomes 



^(OncOj = ^ 
i=l 

where 



fobs{zi) — fth(zi',Q m O>Qr c ) 



E 



H bs{Zi) — H t h(Zi\ flmO, OJ 



H -72 
0.08 



rr'2 tt2 

H th - H o 



((JWl ± zf + OJ 1/2 ± Qli 2 ) 2 + Oo(i ± zf 



(23) 



(24) 



With the same reason as the case of ACDM we take / = 1 as the initial value at high enough redshift. The result is 
illuminated in fig. 2. 

Observing Table 1 carefully, one may find that the datum at z = 3.0 is odd in some degree. From ([T]) we see that in 
ACDM O or ilgie should be smaller than if we require / > 1. Hence our present universe will be curvature dominated 
or becomes an anti-de Sitter (AdS) space, since dust matter redshifts much faster than curvature or vacuum energy. 
Here we give a simple exmple of this problem. In the spatially flat ACDM model, / = 1.46 yields, 



n 6 J n = i.46. 



(25) 



Then we derive Vt m = 2.00, r^e = —1-00. The universe will brake and then start to contract at z = 2.17, which 
completely contradicts to the observations of expansion. So we present fig. 3, which displays the constraint on 
fi m o, Oo m ACDM by H(z) data and growth data, which only include 5 points. Similarly, we plot fig. 4, which 
illuminates the constraint on f2 m o, O c in DGP by H(z) data and growth data, which only includes 5 points. The 
datum at z = 3.0 is excluded. 

Comparing fig. 1 with fig 3., we find the profiles of the two figures are almost the same, but the minimal x 2 , Xmia 
decreases from 12.26 to 9.479. Similarly, comparing fig. 2 with fig. 4, we find the profiles of the two figures are almost 
the same, but Xmin decreases from 12.11 to 9.375. Without the point at z = 3.0 the constraints on Oo 01 ACDM, 
Ono and O c of DGP become more tighten instead. This is also a signal that the datum z — 3.0 may not be well 
consistent with other data. 

If we only consider Xmini we may conclude that DGP is more favored. However, Xmin i s om y °ne point and the 
difference is tiny between the two models. We need more comparisons with the independent results, especially the 
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FIG. 1: 68% , 95% and 99% confidence contour plot of QmO, in ACDM by the growth data in Table I and H{z) data in 
Table II. For 1 a level, Q m0 = 0.275±g;g||| , Q k0 = 0.065±g;itl- 
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FIG. 2: 68% , 95% and 99% confidence contour plot of £l m o, Qko in DGP by the growth data in Table I and H(z) data in 
Table II. For 1 a level, fi m0 = 0.350±o;Jg7 4 , Q, Tc = 0.200lo.o483 ■ 



permitted parameter internals, fitted by the expansion data, which were thoroughly studied. The latest results are 
shown as follows. For ACDM model, ft m0 = 0.279 ± 0.008, Q k0 = -0.0045 ± 0.0065, which are derived from the 
joint analysis of the CMB (five-year WMAP data), the distance measurements from the Type la SN, and the Baryon 
Acoustic Oscillations (BAO) in the distribution of galaxies [H[. For DGP model, fl m0 = 0.28±o°2, ^r c = 0.13 ±0.01 
(SN(new gold)+CMB+SDSS+gas), and Q. m0 = 0.21 ±0.01, 0- c = 0.16 ± 0.01(SN(SNLS)+CMB+SDSS±gas) [H. 

For ACDM, the result of joint fitting by growth data and H{z) data Q m0 = 0.275±^4g, fl k0 = 0.065±^ 4g \ almost 
coincides with the result by expansion data Q m o — 0.279 ± 0.008, flko = —0.0045 ± 0.0065. These types of data are 
well consistent in frame of ACDM model. 

For DGP, the result of joint fitting by growth data and H (z) data impose fl m0 = 0.350±o;£g7 4 , tt rc = O^OOlo ^gg, 
which is not well consistent with the result by expansion data Q m o — 0.281qq2, ^r,. = 0.13 ± 0.01 (SN(new 
gold)+CMB+SDSS±gas), and n m0 = 0.21 ± 0.01, n rc = 0.16 ± 0.01(SN(SNLS)+CMB+SDSS+gas). Concretely 
speaking, Q rc — 0.200^Q'Q4gg(growth-|-iJ(z)) inhabits beyond 2a level of expansion data f2 rc = 0.13 ± 0.01 
(SN(new gold C )±CMB+SDSS±gas). For the data set SN(SNLS)±CMB+SDSS±gas, the result by growth and H{z) 
^mo = 0.350to'Jg74 dwells beyond 3a level of f2 m0 = 0.21 ± 0.01. Therefore DGP model can not fit the observations 
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FIG. 3: 68% , 95% and 99% confidence contour plot of f2 m o, fifco in ACDM by the growth data in Table I and H(z) data in 
Table II. For 1 a level, Q m0 = 0.270±g;^|f , fl k0 = 0.080tg;l£S. The point z = 3.0 in the sample of the growth data is excluded. 
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FIG. 4: 68% , 95% and 99% confidence contour plot of f2 m o, fifco in DGP by the growth data in Table I and H(z) data in Table 
II. For 1 a level, Q m o = 0.345to'Jg4o, ^r c = 0.198ig'.o47i- The point z = 3.0 in the sample of the growth data is excluded. 

of expansion and growth very well at the same time. 

Through the above discussions, we see that the dark energy model is more favored than the DGP model by the 
present data, and the growth data can be an effective probe to study the nature of the dark energy. 

We plot the best fit curves of growth / by growth+H(z) data and expansion data in ACDM, respectively in fig. 
5. Fig. 6 illuminates the best fit curves by growth+H(z) data and expansion data in DGP. It is clear that the gap 
between the best fit curves of growth data and expansion data is much bigger in DGP model than the gap in ACDM 
model. 



IV. CONCLUSIONS 



Perturbation growth is a newly developed method to differentiate between dark energy and modified gravity. In the 
previous works people concentrate on the approximate analytical value of the perturbation growth index of a model, 
and then compare with the observations. But, the index is not a constant in the history of the universe. We fit dark 
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FIG. 5: The perturbation growth / with error bars in Table 1 and the best fit curves in ACDM model. The best fit result by 
growth data and H(z) data inhabits on the red solid curve, the best fit result by growth data except the point z — 3.0 and 
H (z) data resides on the blue dashed one, and and the best fit result by joint analysis of WMAP, SN and BAO dwells on the 
green triangle ones. 
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FIG. 6: The perturbation growth / with error bars in Table 1 and the best fit curves in DGP model. The best fit result by 
growth+H(z) data inhabits on the red solid curve, the best fit result by growth data except the point z = 3.0 and H{z) data 
resides on the blue dashed one, the best fit result by joint analysis of SN(new gold), CMB, SDSS, and gas dwells on the navy 
star curve, and the best fit result by joint analysis of SN(SNLS), CMB, SDSS, and gas is denoted by the pink triangle ones. 



energy and modified gravity models by using the exact evolution equation of perturbation growth. 

The sample of presently available growth data is quite small and the error bars are rather big. Furthermore, we 
always assumed ACDM model for deriving the growth data. Thus it seems proper to fit a model by jointing the 
growth data and the expansion data. The direct H (z) data are new type of data, which can be used to explore the 
fine structures of the Hubble expansion history. We put forward a joint fitting by the growth data and H(z) data. The 
results are summarized as follows: For ACDM, O m0 = 0.275±g;^t ^fco = 0.0651^; For DGP, n mQ = 0.350±g:^|| 4 , 

o — n onn+ ' 0631 

Sir c — U.ZUU_q 0483' 

The minimal % 2 are 12.26 and 12.11 for ACDM and DGP, respectively. The permitted parameters of ACDM by 
growth+H(z) data show an excellent consistency with the previous results inferred from expansion data. However, 
for DGP model the discrepancies of the results of growth+H(z) data and expansion data are at least 2a level. Hence 
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9 in the sense of consistency, ACDM is more favored than DGP. 
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